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This  paper  presents  the  derivation  of  a  linear  model  with  changeable  parameters  of  constant  frequency 
peak  current  mode  controlled  DC/DC  boost  converter  supplied  by  a  PEM  fuel  cell  stack.  The  derived 
model  has  the  same  structure  irrespective  of  the  conduction  mode  and  therefore,  it  is  suitable  for  design 
of  simple  and  advanced  controllers  of  the  output  voltage.  Experimental  results  on  the  system  with  450  W 
boost  converter  supplied  by  PEM  fuel  cell  emulator  show  that  the  derived  model  accurately  describes 
the  system  in  a  given  operating  point,  determined  by  the  load  resistance  or  output  current. 
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1.  Introduction 

Hydrogen  is  a  promising  energy  source  in  the  near  future,  since 
it  can  be  transformed  into  electrical  energy  using  fuel  cells.  All  ben¬ 
efits  of  using  hydrogen  as  possible  replacement  for  fossil  fuels  are 
well  known.  Hydrogen  is  the  cleanest,  the  lightest  and  the  most  effi¬ 
cient  fuel.  Since  hydrogen  is  not  a  primary  energy  source,  it  must 
be  generated  from  other  energy  sources  (e.g.  steam  reformation  of 
natural  gas)  or  by  using  other  energy  sources  for  water  electrolysis. 
Therefore,  other  sources  of  energy  should  be  renewable  (e.g.  wind, 
sun, . . .)  to  remain  in  the  spirit  of  clean  energy  production. 

There  are  numerous  researches  in  which  different  systems 
topologies  with  fuel  cells  are  presented  [  1  -1 1  ].  All  topologies  have 
a  DC/DC  converter  as  the  common  part  of  the  system,  which  sta¬ 
bilizes  the  DC  voltage  from  the  fuel  cell  at  a  certain  (usually  fixed) 
value. 

Control  of  the  fuel  cell  voltage  is  possible  by  controlling  the  air¬ 
flow  through  the  cathode  side  and  the  hydrogen  pressure  on  the 
anode  side  of  the  fuel  cell.  Since  the  fuel  cell  is  a  nonlinear  control 
system,  some  advanced  control  techniques  should  be  used  to  obtain 
the  desired  transient  behavior  of  the  system  and  steady  state  accu¬ 
racy.  In  [12]  an  indirect  model  reference  adaptive  control  is  used 
to  control  the  fuel  cell  voltage,  while  in  [13]  the  model  reference 
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adaptive  control  with  signal  adaptation  algorithm  is  used  for  the 
same  purpose.  A  time  delay  control  is  used  in  [14],  to  improve  the 
fuel  cell  transient  responses  and  at  the  same  time  to  prevent  flood¬ 
ing  and  air  starvation.  Some  advanced  nonlinear  control  techniques 
are  suggested  in  [15-21].  Irrespective  of  the  used  advanced  control 
system,  the  dynamics  of  such  a  system  is  rather  slow  (few  seconds 
or  even  tens  of  seconds),  due  to  the  nature  of  the  control  valves 
and  the  hydrogen  and  airflow  dynamics.  On  the  other  hand,  the 
electrical  part  of  the  fuel  cell  has  a  fast  dynamic  response,  which 
could  not  be  efficiently  controlled  by  input  control  valves.  There¬ 
fore,  it  is  inevitable  to  control  the  DC/DC  converter’s  output  voltage 
to  achieve  fast  transient  responses  (milliseconds)  and  steady  state 
accuracy. 

This  paper  deals  with  the  concept  of  DC/DC  boost  converter  sup¬ 
plied  by  the  fuel  cell.  The  output  voltage  of  the  fuel  cell  BCS  64-32 
[22]  ranges  from  19  V  to  31 V  and  the  boost  converter  has  to  raise 
it  to  50  V. 

The  boost  converter  is  a  nonlinear  and  non-minimum  phase 
system,  and  because  of  that,  its  control  must  be  considered  very 
carefully.  There  are  two  essential  control  techniques  based  on  pulse 
width  modulation  (PWM):  voltage  mode  control  (VMC)  and  current 
mode  control  (CMC).  In  applications  that  require  vary  fast  transient 
responses  the  current  mode  control  is  usually  chosen.  There  are 
three  main  advantages  of  current  mode  control: 

(i)  faster  response  to  the  change  of  referent  voltage  due  to  the 
very  fast  current  loop,  which  acts  like  a  controllable  current 
source, 
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(ii)  the  system  order  is  reduced  by  one,  which  enables  easier 
design  of  the  outer  control  loop  and 

(iii)  inherently  fast  overcurrent  protection. 


The  main  disadvantage  of  the  current  mode  control  is  the  well 
known  instability  for  duty  cycles  greater  than  0.5.  This  problem 
is  easily  solved  by  introducing  a  compensation  ramp  for  constant 
frequency  applications  [23]  or  constant  on/off  times  for  variable 
frequency  applications  [24,25]. 

In  both  control  modes  the  converter  can  operate  in  continuous 
or  discontinuous  conduction  mode. 

Modeling  of  the  VMC  converter  in  continuous  conduction  mode 
is  easily  done  by  using  the  well  known  averaging  techniques 
[26,23].  In  [27]  an  averaged  modeling  of  converters  operating  in 
discontinuous  conduction  modes  is  described. 

Modeling  of  CMC  converters  depends  on  the  chosen  technique 
of  dealing  with  instability  in  duty  cycles  greater  than  0.5.  Literature 
dealing  with  constant  frequency  current  mode  control  [23,28,29] 
usually  concentrates  on  obtaining  high  accuracy  dynamic  models  in 
continuous  or  discontinuous  conduction  modes.  From  the  control 
point  of  view,  this  high  accuracy  for  individual  conduction  modes 
(especially  at  high  frequencies)  is  not  so  important  as  obtaining 
accurate  steady-state  gain  and  dominant  time  constant.  Further¬ 
more,  transfer  functions  that  are  obtained  for  controller  design  are 
not  directly  applicable  because  they  describe  the  behavior  of  the 
output  voltage  with  respect  to  the  change  of  duty  cycle,  rather  than 
the  referent  (peak)  current  with  included  compensation  ramp,  as 
they  should. 

This  paper  focuses  on  obtaining  a  unique  linear  model  with 
changeable  parameters  of  the  CMC  boost  converter,  which  is  valid 
in  continuous  and  discontinuous  conduction  modes  and  is  suitable 
for  the  design  of  a  voltage  control  loop.  The  model  obtained  is  used 
in  [30]  to  design  an  adaptive  control  system,  able  to  achieve  approx¬ 
imately  the  same  dynamic  behavior  of  the  system,  irrespective  of 
the  conduction  mode  and  resistive  load  at  its  output. 

The  model  of  constant  frequency  peak  current  mode  controlled 
DC/DC  boost  converter  is  described  in  Section  2.  The  fuel  cell  model 
is  described  in  Section  3.  The  integrated  process  model  is  described 
in  Section  4  and  the  experimental  identification  of  control  system 
parameters  is  carried  out  in  Section  5.  Finally,  some  conclusions  are 
given  in  Section  6. 


2.  Modeling  of  peak  current  mode  controlled  DC/DC  boost 
converter 

The  model  of  the  converter  applicable  for  the  control  purpose 
must  describe  the  output  voltage  as  a  function  of  the  referent  cur¬ 
rent  and  input  voltage.  The  model  must  have  a  linear  form  in  the 


constant  operating  point  with  parameters  of  the  model  described 
as  a  function  of  the  operating  point  values. 

The  conceptual  scheme  of  the  DC/DC  boost  converter  in  peak 
current  mode  control  is  shown  in  Fig.  1  [23].  Unlike  the  voltage 
mode  control,  the  inductor  current  signal  is  used  for  modulation, 
which  also  has  a  sawtooth  or  triangular  shape.  The  transistor  is 
switched  on  at  the  beginning  of  each  period  by  setting  the  bistable 
(S  input),  and  switched  off  at  the  time  when  the  inductor  current 
reaches  the  referent  current  value  lr  (bistable  R  input). 

The  dynamic  analysis  of  the  peak  current  mode  controlled  boost 
converter  is  based  on  a  method  of  injected  and  absorbed  current 
[23].  According  to  the  method,  the  inductor  and  the  switching  ele¬ 
ments  (comparator,  bistable,  transistor  and  diode)  shown  in  Fig.  1, 
are  replaced  with  a  switching  cell  (Fig.  2).  The  mean  value  of  the 
current  ID i  is  controlled  by  the  referent  signal  Jr,  so  the  switching 
cell  acts  like  a  controllable  current  source  supplying  the  output  RC 
circuit.  The  basic  idea  of  the  injected  and  absorbed  current  method 
is  to  calculate  the  mean  value  of  the  current  /di  ,  which  creates  a 
voltage  drop  in  the  output  RC  circuit,  and  then  calculate  its  total 
differential.  The  procedure  is  applied  for  the  continuous  and  sepa¬ 
rately  for  the  discontinuous  conduction  mode. 


2.1.  Continuous  conduction  mode 


The  inductor  current  in  the  continuous  conduction  mode  never 
falls  to  zero.  In  steady  state  conditions,  assuming  constant  out¬ 
put  voltage  (infinite  output  capacity)  and  ideal  components,  the 
relation  between  input  and  output  voltage  is  given  by  [23]: 


Vout  =  1 

Vin  1  -D’ 

where: 


(1) 


•  Vin  ~  converter  input  voltage  (V), 

•  Vout  -  converter  output  voltage  (V), 

•  D  -  duty  cycle. 


Fig.  2.  The  dynamic  model  of  the  DC/DC  boost  converter  in  peak  current  mode 
control. 
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To  compensate  for  the  well  known  instability  of  the  boost  con¬ 
verter  for  duty  cycle  values  greater  than  0.5,  a  compensation  ramp 
is  used  [23].  The  waveforms  of  the  boost  converter  with  the  com¬ 
pensation  ramp  in  continuous  conduction  mode  are  shown  in  Fig.  3. 

According  to  Fig.  3,  successive  variations  of  current  can  be  deter¬ 
mined  by: 

A/n=f_fcm)n.A/o, 

V  mi  +  my 

where: 


From  (4)  the  duty  cycle  D  follows: 


fr  —  fm 

Vin  +  mL ' 


t  h 


•  m  -  slope  of  the  compensation  ramp  (As-1 ), 

•  mi  -  slope  of  the  inductor  current  rising  part  (As-1 ), 

•  m2  -  slope  of  the  inductor  current  falling  part  (As-1 ), 

•  A/0  -  initial  variation  of  the  inductor  current  (A), 

•  A In  -  nth  successive  variation  of  the  inductor  current  (A), 

•  n  -  ordinal  number  of  successive  variation. 


(5) 


The  current  variation  A In  vanishes  with  time,  i.e.  the  converter 
is  stabilized,  if  the  expression  under  nth  power  is  less  than  1,  or  if: 


m  > 


m2  -  mi 

2 


(3) 


The  inductor  current  k  and  diode  current  im  waveforms  in  the 
continuous  conduction  mode  with  compensation  ramp  are  shown 
in  Fig.  4. 

The  slope  of  the  inductor  current  iL  rising  part,  according  to  Fig.  4, 
is  determined  by: 


mi 


lr  —  f  m  Vfn 

DT  L  ’ 


(4) 


where: 


•  Ir-  referent  inductor  current  (A),  or  control  signal  (mA)  or  (V), 

•  Im  -  minimal  inductor  current  at  the  beginning  of  switching 
period  (A), 

•If  -  minimal  inductor  current  at  the  end  of  switching  period  (A), 

•  T  -  constant  switching  period  (s), 

•  L  -  inductor  inductance  (FI). 


Fig.  4.  The  waveforms  of  the  inductor  current  k  and  diode  current  i'di  in  continuous 
conduction  mode  with  compensation  ramp. 
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The  slope  of  the  inductor  current  iL  falling  part,  according  to 
Fig.  4,  is  determined  by: 

Ir  —  If  inD  —  Vout  —  Vin 


m2  = 


(6) 


(7) 


(l-D)T  3-D  L 

From  (6)  the  expression  for  If  follows: 

If  =  Ir  ~  {(Vout  -  VJnXl  -  D)  -  mDT. 

Inserting  (5)  into  (7)  results  in: 

If  =  Ir  ~  {(Vout  ~  Vin)  +  ^Vout  ~  Vin  -  ml ).  (8) 

The  change  of  minimal  current  value  Jm,  according  to  Fig.  4,  is 
described  by: 

dim  _  If  -  Im 
~dt  ~ 


(9) 


Inserting  (8)  into  (9),  results  in: 

T  =  -'"'-I'"--1 ^ 

The  total  differential  of  change  of  Im  is  given  by: 

d =  k2\  •  dir  +  k22  •  dim  +  k23  *  dVjn  +  /<24  •  dM out ,  (1 1 ) 

where: 


/<21  = 


^22  = 


d(dlm)/(dt)  1  Vout 


dir  T  Vin  “F  PflL 

d(dlm)/(dt)  1  V0ut 


din 


T  Vin  +  ml  ’ 


h  _  d(dlm)/(dt)  _  1  1  Vout  tr  r  ^ 

R 23  —  - tt— -  —  7 x - 7  Ur  —  Im), 

dvin  1  T(Vin  +  mL )2 


(12) 


/<24  = 


d(dlm)/(dt )  _  1  Ir  —  Im  I 


dVout  T  Vm  +  mL  L 

Replacing  the  differential  operator  d  in  (11),  by  small  changes 
A  arround  the  operating  point,  and  applying  Laplace  transform  on 
expression  (11),  results  in: 

s  •  A im(s)  =  /<2i  •  Air(s)  +  /<22  •  Aim(s)  +  k23  •  Ai/in(s) 

+  k24  •  Av0ut{$),  (13) 

from  which  the  expression  for  A im  follows: 

k2\  ■  Air{s)  +  k23  •  A vin{s)  +  /<24  •  A vout(s) 

A'm(S)  “  S^k22  ' 

The  mean  value  of  the  diode  current  ID i  is  given  by  (Fig.  4): 

/D1  =  {  J  b i  (t)dt  =  1(1  -  D)(/r  -  mDT  +  If). 

On  substituting  (5)  and  (8)  into  (15),  the  expression  for  mean 
value  of  the  diode  current  is  obtained,  which  is  also  the  mean  value 
of  the  output  RC  circuit  current: 

j  j  T  L  Ir  —  Im  Vout  -  Vin  -  mL 

/D1  =Ir-  Tl(Vo ut  -  Vin)~  f  Ir-  y— ^  ~  ~L 


(14) 


(15) 


x  (lr  ~  Im)  -  _  Imf 


2  T 


(Vin  +  mLf 


(16) 


by: 


The  total  differential  of  the  mean  diode  current  Jdi  is  determined 


where: 

dhn  L  2Jr-Jm  V0Ut-Vin-mL  L  V0Ut-Vin-2mL 

25  dlr  TVin+mL+  Vin+ml  T  (yin+mLf  " 

,  _  dlm  _L  Ir  V0ut~  Vin  -mL  L  Vout-Vin  -  2mL 

^26  ^  t  \/  ,  ™r  \/  i  mi  t  '  r 


9/m  T  Vin+mL  Vin  +  mL  T  (yirj  +  ml )2 

,,  T  ,  /r-/m  A,  f,  \  ,  L  fI  r  ,2 

*27  —  ~  —  ^77 H  T  I  Tflr—Vout  I  +  —{Ir—lm J 

3V,n  2L  (Vin  +  mL)2  \T  )  2T 

.  d/m  T  Ir-In  L 

^28=7(77 - ^77  +  7 


2Vout-Vin-3mL 

(Vin+mL)3 


(18) 


dV0Ut  21  Vin+mL  2T  (y.n+mQ2 

Replacing  the  differential  operator  d  in  (17)  with  small  changes 
A  arround  the  operating  point,  and  applying  Laplace  transform  to 
(17),  results  in: 


Ai’di(s)  =  k25  •  Air{s)  +  k26  •  A im(s) 

+  /<27  •  Av,-: n(s)  +  k23  •  A V0ut{s). 

Substituting  (14)  into  (19)  yields: 
Ai’di(s)  =  ^/<25  +  •  Alr(s)  +  (j<27  f 

x  •  A Vjn(s)  +  ^/<28  +  5  ^ 


(19) 


^23  ^26  \ 
S  -  /<22  / 


The  change  of  output  voltage  Atw  is  given  by: 
1  +  sRcC 


Avout(s)  =  Aim(s)-R- 


1  +s(Rc  +  R)C‘ 


(20) 


(21) 


On  substituting  (20)  into  (21)  and  rearranging,  the  following 
expression  is  obtained: 


Avout(s)  =  GI)Cl(s)  •  A  iris)  +  Gjc2(s)  •  Avin(s), 
where: 


Gi,cl(s)  = 
Gj,c2(s)  = 

and: 


Avoutis) 


Airis) 

Avoutis) 


AVinis) 


_  k32s 2  +  k33s  +  /<34 
Avin=0  k29S2  +  /<30s  +  /<31  ’ 

_  k33S2  +  /<36S  +  /<37 

Air=0  k2gS2  +  /<30S  +  /<31  ’ 


(22) 


(23) 


^29  =  (i^c  +  R)C  +  k2gRRcC, 

k3o  =  I  —  /<22G(Rc  +  ^)  —  k23R  +  RRcCik22k2g  -  k24k2g), 

/<3i  =  Rik22k28  -  k24k26)  -  k22, 
k32  =  k23RRcC, 

k33  =  RRcG(/<2i^26  -  ^22^25)  +  k23R,  (24) 

L34  =  R(/<2l/<26  -  k22k2s), 

k35  =  k27RRcC, 

k36  =  RRcG(/<23^26  -  k22k27)  +  k27R, 
k37  =  Ri k23 k2Q  -  k22k27). 

By  neglecting  the  equivalent  series  resistance  of  the  capacitor 
Rc ,  the  transfer  functions  (23)  assume  the  form: 


Gj,cl(s)  —  I^i,c\  ' 

Gj,c2(^)  —  I^i, c2  ' 

where: 
v  j<34 

I<Lcl  =  k7/’ 


I  Ti.cb/S 


1  +  2tiTi  cS  +  T?cs2  ’ 

I  T  Ti.cb2$ 

1  +  2?,Tj,cS  +  T^s2  ’ 
17  fc37 

K'c2  ~  hi' 


(25) 


x  fc33  T  fc36 

Ti’cM_  k34’  Ti'cb2~k31' 

J  /^29  r  _  ^30 

!’c  _  V  k31  ’  1  “  2V& 


(26) 


31 


d/j)i  =  l<25  ■  dir  +  ^26  •  d/m  +  ^27  '  dVln  +  k2 8  •  dVout , 


(17) 
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By  the  choice  of  a  much  higher  switching  frequency  than  the 
frequency  of  the  output  RC  circuit,  the  transfer  function  (25)  poles 
become  real,  i.e.  &  >  1,  so  the  transfer  functions  (25)  can  be  written 
as: 


Q,ci  (s)  —  I<i,cl 
Q,c2(s)  =  fQ,c2 


1 

(1  +  7f,cls)(l  +  ^i,c2s) 

l+hcb2l 

(1  +  Tf,cls)(l  +  Tj  C2s) 


(27) 


where: 


he  1  = 
7;,c2 


f 

'  i,  c 

&  +  \/^!?  “  1 


(28) 


Furthermore,  if  >  1  then  from  (28)  it  follows  that  TIiCl  »  T2-c2, 
so  Ti  c2  can  be  neglected  with  respect  to  TI  cl.  The  nonminimum- 
phase  zero  in  transfer  function  G2iCi  (27)  can  be  neglected  only  if  a 
relatively  small  value  of  the  inductor  inductance  L  is  chosen.  The 
zero  in  transfer  function  Gic2  (27)  can  only  be  neglected  if  the  equiv¬ 
alent  series  resistance  of  the  output  capacitor  is  neglected  (ceramic 
capacitor  in  parallel  with  electrolyte  capacitor). 

It  follows  that  by  the  choice  of  the  inductance  of  a  relatively 
small  amount,  the  transfer  functions  (27)  are  simplified  to  the  PTi 
form: 


Q,ci  (s) 

Q,c2(s) 


Kj,cl 

1  +  TUaS 

I<i,c2 

1  +  7),  cl 5 


(29) 


2.2.  Discontinuous  conduction  mode 


The  inductor  current  falls  to  zero  during  one  switching  cycle  in 
the  discontinuous  conduction  mode.  The  relation  between  input 
and  output  voltage  in  steady  state  conditions,  assuming  constant 
input  and  output  voltages  and  ideal  components,  is  given  by  [23]: 


Vout  _  Di  +  D 

v~~~d^~ 


(30) 


where:  D i  -  a  fraction  of  the  switching  cycle  in  which  the  inductor 
current  is  zero. 

The  waveforms  of  inductor  current  iL  and  diode  current  im  in 
the  discontinuous  conduction  mode  with  the  compensation  ramp 
are  shown  in  Fig.  5. 

The  slope  of  the  inductor  current  iL  rising  part,  according  to  Fig.  5, 
is  determined  by: 


mt  =  -m  = 

1  DT  L  ’ 


from  which  the  expression  for  D  follows: 


(31) 


Fig.  5.  The  waveforms  of  inductor  current  k  and  diode  current  io i  in  discontinuous 
conduction  mode  with  compensation  ramp. 


The  diode  current  mean  value  ID i  is  determined  by  (Fig.  5): 

/D1  =  ^(fr-mDT).  (36) 

Substituting  (32)  and  (35)  into  (36),  results  in: 


A  1 

2  T  Vout-Vin 


Vl 


(Vin  +  mL)2' 


(37) 


The  total  differential  of  the  diode  current  mean  value  ID i  has  the 
form: 


dlo  1  =  /<38  dir  +  l<39dVin  +  k4odV0ut, 


(38) 


d=4- 


h 


(32) 


T  Vin  +  mL ' 

The  slope  of  the  inductor  current  ii  falling  part,  according  to 
Fig.  5,  is  determined  by: 

Vout  -  vin 


lr  D 

m2  =  7^7  -  m  •  — 
1  D{Y  Di 


(33) 


where  D\  is  marked  in  Fig.  5.  From  (33)  the  expression  for  D \  fol¬ 
lows: 


L  lr-  mDT 
A~T'V0Ut-Vin 

Substituting  (32)  into  (34),  results  in: 
lr  Vin 


(34) 


T  Vout  —  Vfn  Vin  +  mL 


where: 

,  dJm  L  vl± 

Mr  T  (Vout  -  Vin)(Vin  +  mL)2  ’ 
dip ,  LVmlj  Vfn  +  rnL(2V0ut  -  V,„) 

dvin  2 T  (Vout  -  y,'n)2(Vm  +  mL)3  ’ 

c)fD1  L 

3 Vout  2T  (V!rl  +  mL)2(V0U[  -  Vjn)2  ' 


(39) 


Replacing  the  differential  operator  d  in  (38)  with  small  changes 
A  arround  the  operating  point,  and  applying  Laplace  transform  to 
(38),  results  in  expression  for  the  small  changes  of  output  RC  circuit 
current  mean  value: 


(35) 


Ai’di(s)  =  /<38  Air(s)  +  /C39  Apin(s)  +  1<4oAv0ut{s). 


(40) 
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The  change  of  output  voltage  A vout  is  obtained  by  multiplying 
AiDi  and  output  circuit  impedance: 

(4,) 

On  expressing  the  AiD1  from  (41 )  and  inserting  it  into  (40),  the 
change  of  output  voltage  is  obtained: 


A Vout(s)  —  Q,dl(s)  •  Al’r(s)  +  Gi  d2{s )  •  A Vjn(s), 


where: 
Q,dl(s)  — 
Q,d2(s)  = 


A  Vout(s) 
Air(s) 
A  t'ont(s) 
AVin(s) 


I<Ucn 

1  +  ^i,ds 
*i,d2 
1  +  ^i,ds 


Ayin=0 

Air=0 


(42) 


(43) 


are  the  appropriate  transfer  functions  in  the  discontinuous  conduc¬ 
tion  mode,  with  parameters: 

^  k38R(\+RcC) 

Ki'dl  =  1  -k40R  ’ 

„  k39K(l+KcC) 

Ki  d2  =  1  -  k40R  ’ 

T  ftc+ft  ^ 

l-k40R 


(44) 


It  is  important  to  notice  that  with  careful  converter  design  the 
transfer  functions  for  continuous  (29)  and  discontinuous  (43)  con¬ 
duction  modes  have  the  same  PTi  form,  with  different  parameters. 


•  I  pc  -  fuel  cell  current  (A), 

•  §i_4  -  cell  parameters,  whose  values  are  defined  based  on 
theoretical  equations  with  kinetic,  thermodynamic  and  electro¬ 
chemical  foundations  [33], 

•  Co2  -  oxygen  concentration  in  the  catalytic  interface  of  the  cath¬ 
ode  (mol cm-3),  defined  by: 


5.08  106  -  exp(-498/T)’ 


The  parameter  £2  is  determined  by  [31  ]: 
i 2  =  0.00286  +  0.0002  ln(A)  +  4.3  •  10"5  ■  ln(CH2 ),  (52) 

where  Ch2  is  determined  by: 


Ch2  = 


_ PH2 _ 

5.08  •  10s  •  exp(-498/T)' 


(53) 


The  resistive  or  the  ohmic  losses  result  from  the  resistance  to 
proton  transfer  through  the  membrane  and  the  resistance  to  the 
electron  transfer  through  the  electrodes: 


Vohm  =  he  ■  {R-M  +  Rc)  =  he  ■  Rohm >  (54) 

where  Rc  is  the  resistance  to  the  electron  trransfer  through  the 
electrodes,  usually  constant,  and  RM  is  the  membrane  equivalent 
resistance,  defined  by  [31  ]: 

(55) 

where: 


3.  Modeling  of  PEM  fuel  cells 

The  fuel  cell  is  an  electrochemical  device  that  directly  converts 
the  chemical  energy  of  the  fuel  (usually  hydrogen)  into  electric¬ 
ity.  The  byproducts  of  the  reaction  are  water  and  heat.  The  basic 
electrochemical  reactions  are  described  by: 

•  Anode: 


•  pm  -  specific  resistance  of  the  membrane  (^  cm), 

•  A-  cell  active  area  (cm2), 

•  l  -  membrane  thickness  (cm). 


Membranes  of  the  Nafion  type,  considered  here,  are  registered 
trademark  of  the  Dupont  company,  and  are  widely  used  in  the  PEM 
fuel  cells.  Their  specific  resistance  is  given  by  [33]: 


181.6  [1  +  0.03  •  Ifc/A  +  0.062  ■  (T/303)2{IFC/A)2-5] 
[Xm  -  0.634  -  3  ■  IFC /A\  ■  exp[4.1 8  •  (T  -  303/T)] 


H2  ->■  2H+  +  2e~,  (45) 

•  Cathode: 

lo2  +  2H+  +  2e-  H20,  (46) 

•  Overall  reaction: 

H2  +  lo2  -  H20.  (47) 

The  single  cell  voltage  can  be  described  by  [31  ]: 

Vcell  =  E-Vact-Vohm-Vcon.  (48) 

The  thermodynamic  or  the  Nernst  potential  E  of  the  cell  repre¬ 
sents  its  reversible  voltage,  and  is  given  by  [32]: 

E  =  1 .482  —  8.45  ■  10-4  •  T  +  4.31  ■  10-5  •  T  ■  ln(pH2  (49) 

The  voltage  drop  due  to  the  activation  of  the  anode  and  cathode, 
also  known  as  the  activation  overpotential,  is  given  by: 

Vact  =  -l^+&-T  +  $3-T.  ln(Co2 )  +  §4  •  T ■  ln(/FC)],  (50) 


where  Xm  is  the  parameter  which  determines  the  average  water 
content  in  the  membrane,  with  reported  values  arround  14-23  in 
conditions  of  100%  membrane  humidity. 

The  concentration  polarization  or  the  voltage  drop  due  to  mass 
transport  occurs  when  the  reactant  is  rapidly  consumed  in  the 
electrochemical  reaction,  so  it  is  necessary  to  quickly  deliver  more 
reactant  to  the  reaction  site.  Given  the  limited  permeability  of  the 
electrodes,  there  is  a  concentration  gradient  of  the  reactant.  This 
occurs  at  higher  currents,  because  it  is  known  that  the  amount  of 
hydrogen  consumed  (in  mols  per  second)  is  directly  proportional 
to  the  fuel  cell  current: 

"h2  =  (57) 

Therefore,  the  maximum  current  density  Jmax  (A cm-2)  is 
defined,  which  corresponds  to  the  maximum  delivery  rate  of  the 
reactant  in  the  reaction  (catalytic)  layer.  The  fuel  cell  is  designed 
so  as  to  have  as  great  a  value  of  maximum  current  density  Jmax  as 
possible.  Typical  values  range  from  0.5  to  2  A  cm-2,  while  carefully 
designed  fuel  cells  can  provide  up  to  3  A  cm-2. 

Thus,  the  voltage  drop  due  to  concentration  polarization  can  be 
associated  with  the  current  density  in  the  following  way  [31  ]: 

Vcon  =  — B  •  In  (l  —  t^— )  =  — B  •  In  (l  —  >  (58) 

V  Jmax  J  \  Jmax'7\J 


where: 


where: 
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Fig.  6.  The  equivalent  electrical  scheme  of  the  dynamic  fuel  cell  model. 

•  B  -  fuel  cell  parameter  (V), 

•  J  -  fuel  cell  current  density  (A cm-2). 

The  fuel  cell  dynamics  can  be  described  by  first  order  delay  in 
the  activation  and  concentration  polarization  [31  ],  with  associated 
time  constant: 

T  =  CRa,  (59) 


Since  changes  in  pressure  of  the  reactants  in  the  models  of  volt¬ 
age  losses  are  taken  into  account  through  semiempirical  relations, 
the  pressures  of  the  reactants  actually  become  the  constants  of  non¬ 
linear  and  linear  models.  Furthermore,  the  voltage  dependence  on 
temperature  is  significant,  but  with  a  dynamic  point  of  view  it  is 
irrelevant,  because  the  temperature  of  the  fuel  cell  must  be  regu¬ 
lated  [36,37],  and  hence  it  can  be  considered  constant.  It  is  therefore 
necessary  to  find  a  linear  model  of  voltage  dependence  on  fuel  cell 
current: 


Vfc  =/(/fc).  (65) 

The  Nernst  voltage  does  not  depend  on  the  fuel  cell  current,  so 
it  can  be  written  as: 

E=gi  =  1.482  -  8.45  10“4  r +  4.31  10“5  r  -  ln(pH2  pg25).  (66) 

The  activation  voltage  drop  (50)  depends  on  the  fuel  cell  current, 
so  the  expression  (50)  can  be  written  as: 

Vact  =  g2+ S3  In  {Ifc),  (67) 

where: 

52  =  ~(§i  +  &T  +  §3Tln(Co2)), 

53  =  -§4  T, 

where  parameter  £2  is  determined  by  expressions  (52)  and  (53), 
and  concentration  C02  is  determined  by  (51). 

The  concentration  polarization  also  depends  on  the  fuel  cell  cur¬ 
rent  (58),  so  it  can  be  written  as: 

Vcon=g4ln(l+g5/Fc),  (69) 

where: 


54  =  ~B, 

55  =  ~ 


1 


Jmax  •  A 


(70) 


where  C  represents  the  cell  equivalent  capacitance  in  Farads,  while 
Ra  represents  the  equivalent  resistance  in  £2.  The  capacitance  value 
C  usually  has  a  constant  value  of  several  Farads,  while  the  resistance 
value  can  be  determined  as: 


From  (54),  (55)  and  (56),  ohmic  voltage  drop  can  be  written  as: 


Vohm  =  he 


gl+g8lFC+g9ljl\ 

gw+gnhc  J 


(71) 


Ra  =  Vact,+  V"con.  (60) 

he 

Fig.  6  shows  the  equivalent  electrical  scheme  of  the  dynamic 
fuel  cell  model.  From  Fig.  6  it  follows: 


V cell  ~  B  Vc  Vohm  > 


(61) 


where  Vc  is  the  capacitor’s  C  voltage  (Fig.  6),  determined  by  the 
differential  equation: 


dVc  _  he  _ Vc  \ 

dt  C  \  Vact  +  VCon  ) 


(62) 


The  voltage  of  the  whole  stack  is  determined  by: 


VFC  =  Nc.Vcelh 


(63) 


where  Nc  is  the  total  number  of  cells  in  the  fuel  cell  stack. 
3.1.  Linearized  fuel  cell  model 


where: 
S6  =  Re, 


l 


g7  =  181.6- 

Al 

g8  =  5.448--, 

S3  =  1-2264- 10-4-  JLj, 
gio  =  (Am -0.634).  exp 
3 


Sn  =  -j  ■  exp 


4.18- 


4.18- 

T  —  303  \ 

T  ) 


'T-303' 


Substituting  (67)  and  (69)  into  (62),  results  in: 


dVc 

dt 


he 

C 


1  - 


(72) 


(73) 


g2  +  S3 1  n(hc )  +  S4  ln(  1  +  Sshc )  > 

The  change  of  voltage  Vc  (73)  is  a  function  of  two  variables,  IFC 
and  Vc  itself,  so  its  total  differential  is: 


The  linearized  fuel  cell  model,  suitable  for  controller  design 
purposes,  is  derived  and/or  described  in  [1,34,35].  Due  to  the  com¬ 
pleteness  and  some  minor  differences,  the  linear  model  will  be 
derived  once  again  in  the  following  lines. 

From  the  previously  derived  equations,  it  is  evident  that  the 
voltage  of  the  fuel  cell  is  a  function  of  several  variables: 

VFc=nhe,T,pH2,p02).  (64) 


d~lf  =gi2  -dlfc  +gi3  dVc,  (74) 

where: 


d{dVc)/{dt)  _  Ipc_  (S3 )/ (he )  +  (S4S5 )/( 1  +  Sshc ) 

3  he  B  g2  +  S3  In  (he )  +  S4  ln(  1  +  Ss  he )  ’ 

d(dVc)/(dt)  _  f _ Ipc _ 

dVc  C  S2+ S3  ln(/Fc )  +  S4  ln(  1  +  Sshc ) 
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Substituting  (66)  and  (71 )  into  (61 ),  and  then  in  (63),  results  in: 

§7  +  &8 he  +  g9lpc  " 


Vfc  =  Nc  gi  -  Vc  -  Ifc  S6  + 


&io  +guhc 


(76) 


The  fuel  cell  voltage  VFC  (76)  is  a  function  of  two  variables,  IFC 
and  Vc,  so  its  total  differential  is: 


dVFc  =  gi4  •  dIFc  +gi5  •  dVc, 
where: 

dVFC  AT  (  §7  H~  §8^fc  +g9ljc5 

gl4=^  =  ~Nc  g 6  +  a  1  ~ 

0lFC  \  gw+guhe 

(gshc  +  2.5ggrir5)(gio  +guhc) 


(77) 


-Ndpc 


+NclFCgU 


(&io  +gnhc) 
gl  T  g8^FC  T  g9^pc 


(78) 


gl5  = 


dVpc 

dVr 


(&io  +gnhc) 
=  —Nr. 


Replacing  the  differential  operator  d  with  small  changes  A 
arround  the  operating  point,  and  applying  Laplace  transform  to  (74) 
and  (77)  and  rearranging,  results  in  the  transfer  function  of  the  fuel 
cell: 


CFc(s)=4^=KFC.1+7k’bS 


A  ifc(s) 


1  +  /fc+s’ 


(79) 


where: 


I<fc  = 


&13&14  -&12&15 


gl3„ 


TFC,b  =  - 


£14 


&13&14  -&12&15 


(80) 


TFC ,1  =  -  —  • 
g  13 


1 

+  2 
1  Vout 

+  2  '  Vin  +  ml  ' 


1  l  (  Ir-Irr 


T  Vin  +  mL\  Vin  +  mL 


(Vout  -  vin  - 


mL)) 


04 


Ir-lm  \ 
Vin  +mLJ 


(82) 


The  mean  value  of  boost  converter  inductor  current  in  the  dis 
continuous  conduction  mode,  according  to  Fig.  5,  is: 

I? 


1  L  VoutVin 
L’d  2  1  Vout  -  Vm  ■  [v.n  +  mLf  ■ 

Parameter  p  for  discontinuous  conduction  mode  {p  =  pd)  is 


Pd 


dh,d  _  L  V0UtVin 


Ir 


(83) 


(84) 


9 Ir  T  Vout  ~  Vin  (Vin  +  mL)2 

The  transfer  function  of  the  whole  system,  according  to  Fig.  7 


is: 


G(5)  =  ^“(s) 3  =  G,'-l(s)  +  p '  Gfc(s)  ‘  C|'2(S)- 


(85) 


The  transfer  function  of  the  fuel  cell  has  the  form  (79),  while  the 
boost  converter  transfer  functions  have  the  PTi  form  (29)  or  (43), 
depending  on  the  conduction  mode: 


G^(S)  =1  +  7^- 


(86) 


Substituting  (79)  and  (86)  into  (85),  results  in  the  transfer  func¬ 
tion  of  the  integrated  system  in  the  operating  point: 


1  +  TjjS 


_  ^^out(s)  _ 

Airis)  '  (1  +1,5X1  +I2S)’ 

where  parameters  in  the  operating  point  are: 


(87) 


4.  Integrated  process  model 


Integration  of  nonlinear  fuel  cell  model  and  the  boost  converter 
is  done  simply  by  connecting  the  inductor  current  signal  lF  to  the 
current  input  of  the  fuel  cell  IFC,  and  connecting  the  fuel  cell  voltage 
signal  VpC  to  the  input  of  the  boost  converter  Vin. 

Unlike  the  integration  of  nonlinear  models,  the  integration  of  the 
linearized  models  necessary  for  designing  the  control  circuit  cannot 
be  implemented  directly.  The  reason  is  that  the  linearized  model 
of  the  boost  converter  does  not  contain  the  inductor  current  +,  but 
the  referent  current  Jr.  It  is  therefore  necessary  to  find  a  connection 
between  the  inductor  mean  current  lF  and  the  referent  current  /r, 
and  include  it  in  the  linearized  model  with  parameter  p  (Fig.  7). 

The  mean  value  of  boost  converter  inductor  current  in  the  con¬ 
tinuous  conduction  mode,  according  to  Fig.  4,  is: 


,  =  1  L  k-lm 

Uc  2  '  T  '  Vin  +  mL 


Ir-h 

2  V  T 

n 


Vm  +  mL 
Ir-h 


i  1  ,  1  ,  Ir  —  h 71 

■Im+  ~h  -  ~mL  -  — - 

2  2  V,n  +  mL 


)  ■  (jr  ~  I(Vout  ~  Vin )) 


)  -*■')■  (8,) 


T  Vin  +  mL  j 

Parameter  p  for  continuous  conduction  mode  (p  =  pc)  is: 


_  9 IF,c  _  I  f  *  L  Im  —  mL  \ 

pc~  w  -  2  V  f  ■  Vin  +  mLj 


l  r  l  L 
+  2  [T  Vin  +  mL 


T 

L 


(.Vout  -  Vin)) 


K  =  l<u  +pKFCKi2, 

J  _  Ki,  1  TfC,  1  +  PTcT  2  Tpc,  b 

b~  I<u  + 2  ’  (88) 

h  =  Tpc,u 
T2  =  Tj. 

After  obtaining  the  dynamic  model  of  the  whole  system,  it  would 
be  convenient  to  run  some  simulations  that  would  confirm  the 
obtained  theoretical  results.  However,  the  simulations  of  the  non¬ 
linear  model  of  the  converter  in  Matlab  -  Simulink  require  very 
small  step  size  in  a  numerical  integration  method.  According  to 
the  authors’  experience,  a  variable  step  solver  should  be  used,  with 
the  maximum  step  size  of  1  If 1200  for  the  discontinuous  and  even 
smaller  for  the  continuous  conduction  mode  (I If 1500),  where /is 
the  switching  frequency  of  the  converter.  This  leads  to  very  slow 
simulations  and  rapid  consumption  of  the  computer  memory.  It  is 
often  impossible  to  get  the  converter  to  the  operating  point  (start¬ 
ing  from  zero  initial  conditions)  and  then  apply  the  request  for  a 
step  change  around  the  operating  point  to  record  the  transient. 
It  helps  to  use  a  specialized  software  (Spice  engine  based),  but 
they  lack  of  support  for  other  components,  like  the  fuel  cell  in  our 
example.  Calculation  of  necessary  initial  conditions  to  be  in  the 
operating  point  immediately  is  possible,  but  only  for  mean  values 
of  the  inductor  current  and  input  and  output  voltages,  so  it  is  still 
inevitable  to  adjust  the  initial  conditions  manually,  which  is  very 
time  consuming  due  to  the  very  slow  simulations.  Therefore,  the 
authors  decided  to  build  the  special  version  of  the  converter  which 
enables  all  the  experiments  that  would  usually  be  done  through 
computer  simulations.  The  experimental  system  and  identification 
of  its  parameters  are  described  in  the  next  section  of  the  paper. 
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Fig.  7.  Block  diagram  of  the  entire  linearized  system. 


Fig.  8.  The  part  of  the  hydrogen  facility  in  the  Laboratory  for  Renewable  Energy 
Sources  -  LARES. 


5.  Experimental  identification  of  the  control  system 
parameters 

The  control  system,  installed  in  the  Laboratory  for  Renewable 
Energy  Sources  (LARES)  [38],  consists  of  a  500  W  fuel  cell  BCS  64- 
32  (Fig.  8),  manufactured  by  BCS  Fuel  Cells  Inc.,  a  450  W  boost 
converter  specifically  designed  and  built  in  the  Croatian  company 
Mareton  Power  Electronics  (Fig.  9),  a  Magna  Power  Electronics  fuel 
cell  emulator  (Fig.  10)  [39,40],  and  finally  a  Compact  RIO  (cRIO) 
9024  digital  controller  with  input-output  modules,  manufactured 
by  the  National  Instruments  company,  which  will  be  used  for  the 
experimental  identification  of  the  system  parameters  and  imple¬ 
mentation  of  digital  control  algorithms  [41,42]. 


Fig.  9.  DC/DC  boost  converter,  designed  in  Mareton  Power  Electronics. 


5.1.  Experimental  identification  of  the  boost  converter  parameters 

A  conceptual  scheme  of  the  boost  converter,  built  in  Mareton 
Power  Electronics,  is  shown  in  Fig.  11. 

The  converter  is  designed  for  the  input  voltage  range  19-31 V 
and  nominal  output  voltage  50  V.  The  switching  frequency  is 
constant  and  set  to  100  kHz.  The  declared  efficiency  is  97%  mini¬ 
mum.  The  nominal  output  current  equals  9  A,  while  the  maximum 
allowed  current  equals  13  A. 

Experimental  identification  of  the  boost  converter  is  carried  out 
in  several  operating  points,  determined  by  the  output  current  Iout- 
The  output  current  of  the  converter  in  the  normal  operating  condi¬ 
tions  ranges  in  the  interval  Iout  e  [W, min.  lout, maxi 

The  first  operating  point  is  determined  by  the  converter  output 
current  Iout  =  I out, min ,  and  the  last  by  Iout  =  lout, max •  The  first  oper¬ 
ating  point  represents  the  minimum  energy  consumption  of  the 
system.  It  is  assumed  that  a  viable  system  design  can  provide  the 
amount  of  power  required  to  supply  additional  components  (com¬ 
pressor  or  blower,  the  cooling  unit  and  a  system  for  humidifying 
gases  if  present)  to  be  about  10  percent  of  the  nominal  power. 
For  the  fuel  cell  model  BCS  64-32  this  is  a  value  of  2.5  A,  where 
the  boost  converter  output  current  is  approximately  Iout,min  =  1  A. 
The  last  operating  point  is  determined  by  the  boost  converter  rated 
current,  or  nominal  load  of  the  fuel  cell  I0ut,max  =  9  A. 

The  output  voltage  is  reduced  to  0-1 0  V  range  by  using  the  resis¬ 
tor  divider  with  coefficient  /<^  =  0.176,  and  the  voltage  feedback 
signal  is  additionally  digitally  filtered  with  time  constant  =  350 
[xs,  to  reduce  the  voltage  ripple  about  10  times  (20  dB).  Therefore, 
the  feedback  transfer  function  has  the  form: 

Avjh(s)  _  Kfl,  _  0.176 

Avout(s)  1+T^s  1  + 0.00035s'  1 1 


Fig.  10.  Front  panels  of  the  fuel  cell  emulator  Magna  Power  Electronics  and  the 
programmable  electronic  load  HP  6050A. 
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NI  9265  cRIO-9024  NI  92 1 5 

Fig.  11.  A  conceptual  scheme  of  the  boost  converter,  built  in  Mareton  Power  Electronics. 


The  transient  responses  for  the  two  described  (boundary)  oper¬ 
ating  points  are  shown  in  Figs.  12  and  13.  The  responses  are 
recorded  in  LabviEW  and  transfered  into  Matlab  for  easier  presen¬ 
tation  and  processing  of  results. 

Changes  around  the  operating  points  are  described  by  the  lin¬ 
earized  model  (86),  while  the  parameters  of  transfer  functions  were 


determined  using  the  simplex  method  of  optimization  in  Matlab, 
with  minimization  of  the  ISE  criterion  error: 

r* 2 

Icrit  —  /  {AVout  ~  Av, out,lin )  (90) 


0  0.2  0.4  0.6  0.8  1  1.2  1.4  1.6 

t( S) 


Fig.  12.  Transient  responses  of  the  output  voltage  Vout  and  feedback  voltage  V on  the  change  of  control  signal  /r(t),  with  Iout  =  1  A,  V/n  =  28  V. 
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Fig.  13.  Transient  responses  of  the  output  voltage  Vout  and  feedback  voltage  on  the  change  of  control  signal  Ir(t ),  with  Iout  =  9  A,  Vin  =  20.5  V. 


where  t\  and  t2  are  initial  and  final  simulation  times,  i.e.  the 
extracted  part  of  the  recorded  transients. 

The  recorded  responses  and  converter  linearized  model  tran¬ 
sients,  along  with  their  deviations,  are  given  in  Figs.  14  and  15. 

After  completing  the  optimization  procedure,  the  transfer  func¬ 
tion  Gpi  for  the  first  (boundary)  operating  point  determined  by 


lout  =  1  A  is  obtained,  as  well  as  the  transfer  function  Gp9  for  the 
last  (boundary)  operating  point  determined  by  Iout  =  9  A: 

,  .  _  Avout(s)  _  12489 

pU  '  Air(s)  1  + 0.0745s’ 

^  _  AiWs)  _  5037.5 

Gp9(S;  “  Air(s)  “  1  +  0.0171s 


s> 

<1 


Fig.  14.  The  recorded  converter  output  voltage  transient  without  initial  value 
and  linearized  model  Gpi  (91)  transient  on  the  change  of  control  sig¬ 
nal  /r(t)  =  5.823  +  0.1S(t- 0.092) -0.1S(t- 0.492) mA,  with  conditions  described  in 
Fig.  12  caption. 


Fig.  15.  The  recorded  converter  output  voltage  transient  without  initial  value 
and  linearized  model  Gpg  (91)  transient  on  the  change  of  control  sig¬ 
nal  /r(t)  =  17.27  +  0.1S(t-0.075)-0.1S(t-0.175)mA,  with  conditions  described  in 
Fig.  13  caption. 
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Table  1 

The  parameters  of  the  linear  model  of  the  boost  converter  (86)  /<,  i  and  T,  obtained 
in  several  operating  points  determined  by  the  output  current  Iout- 


lout  (A) 

Vin  (V) 

% 

Ti  (ms) 

1 

28.0 

12,489 

74.5 

3 

24.7 

7116.5 

22.7 

5 

23.5 

6496.8 

19.5 

7 

22.2 

6022.8 

18.3 

9 

20.5 

5037.5 

17.1 

Fig.  16.  Static  V-I  and  P-I  curve  of  the  fuel  cell  BCS  64-32,  obtained  from  the  man¬ 
ufacturer,  BCS  Fuel  Cells,  Inc.  [22]. 

Figs.  1 4  and  1 5  show  that  the  maximum  error  of  the  linear  model 
is  actually  determined  by  a  standard  deviation  of  the  output  volt¬ 
age  ripple  and  noise,  and  that  the  derived  transfer  functions  (86) 
accurately  describe  the  behavior  of  the  peak  current  mode  DC/DC 
boost  converter. 

The  obtained  results  for  all  operating  points  are  summarized  in 
Table  1. 


The  transient  responses  in  other  operating  points,  given  in 
Table  1,  are  similar  to  the  responses  in  the  described  boundary 
operating  points,  which  means  that  maximum  error  of  the  linear 
model  is  also  determined  by  the  standard  deviation  of  the  out¬ 
put  voltage  ripple  and  noise.  In  explanation,  the  transfer  functions 
(86),  along  with  the  appropriate  parameters  from  Table  1,  accu¬ 
rately  describe  the  behavior  of  the  peak  current  mode  DC/DC  boost 
converter  in  the  given  operating  point  determined  by  the  output 
current. 

5.2.  Experimental  identification  of  the  boost  converter  supplied 
with  the  fuel  cell  emulator 

The  boost  converter  supplied  by  the  fuel  cell  emulator,  with 
stored  static  V-I  curve  according  to  manufacturer’s  data  (Fig.  16), 
is  identified  in  the  same  operating  points  as  the  boost  converter. 
Thereby,  the  fact  that  the  boost  converter  is  already  identified  is 
used,  and  its  time  constants  are  given  in  transfer  functions  (91) 
and  in  Table  1. 

The  responses  of  the  boost  converter  supplied  by  the  fuel  cell 
emulator  for  the  two  described  boundary  operating  points  are 
shown  in  Figs.  17  and  18. 

Changes  around  the  operating  point  are  described  by  the 
linearized  model  (87),  and  transfer  function  parameters  are 
determined  using  simplex  optimization  method  in  Matlab,  by  min¬ 
imizing  integral  square  error  (ISE)  criterion: 

ft 2 

firit  —  /  ~  ^Vfb,lin)  dt,  (92) 

Jt  i 

where  t\  and  t2  are  initial  and  final  simulation  times,  i.e.  the 
extracted  part  of  the  recorded  transients. 

The  recorded  transients  and  responses  of  the  linearized  model, 
together  with  their  deviations,  are  shown  in  Figs.  19  and  20. 

After  optimization  is  carried  out,  the  transfer  functions  are 
obtained,  G\  for  the  first  (boundary)  operating  point  determined  by 


0  0.2  0.4  0.6  0.8  1  1.2  1.4  1.6 

t  (s) 


Fig.  17.  Transient  responses  of  the  boost  converter  supplied  by  the  fuel  cell  emulator  (output  voltage  Vout  and  feedback  voltage  V}&)  on  the  change  of  control  signal  Ir(t ),  with 
W  =  1  A,  Vin  =  VFC  =  28V. 
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Fig.  18.  Transient  responses  of  the  boost  converter  supplied  by  the  fuel  cell  emulator  (output  voltage  Vout  and  feedback  voltage  V^)  on  the  change  of  control  signal  Ir(t),  with 
W  =  9A,V!n  =  VfC  =  20.5V. 


lout  =  1  A,  and  G9  for  the  last  (boundary)  operating  point  determined 
by  lout  =  9  A,  with  included  feedback  element  (89): 

_A  Vfh(s)_  2088  (1  +  0.0186s) 

llSj  “  Air{s)  “(1  +0.0745s)(l  +0.016s)(l  +  0.00035s)’ 

(93) 

=  Avfb(s)=  461.(1  + 0.1576s) 

9lSj  “  Air(s)  “(1  +0.1142s)(l  +0.0171s)(l  +  0.00035s) ‘ 


Figs.  19  and  20  show  that  the  maximum  errors  of  linearized 
models  are  dominantly  determined  by  the  standard  deviation  of  the 
output  voltage  ripple  and  noise,  so  the  derived  transfer  functions 
(87)  accurately  describe  the  behavior  of  the  DC/DC  boost  converter 
supplied  by  the  fuel  cell  emulator. 

The  obtained  results  for  all  operating  points  are  summarized  in 
Table  2.  Additionally,  the  Bode  plot  of  the  frequency  characteristics 
of  the  obtained  model  (87),  with  parameters  given  in  Table  2,  are 
shown  in  Fig.  21. 


Fig.  19.  The  recorded  transient  of  converter  output  voltage  without  initial  value 
and  the  transient  response  of  the  model  Gi  (93)  on  the  change  of  control  signal 
Jr(t)  =  5.823  +  0.1S(t-0.37)-0.1S(t-0.97)mA,  with  conditions  according  to  Fig.  17 
caption. 


Fig.  20.  The  recorded  transient  of  converter  output  voltage  without  initial  value 
and  the  transient  response  of  the  model  Gg  (93)  on  the  change  of  control  signal 
Jr(t)  =  17.27  +  0.1S(t-0.37)-0.1S(t-0.97)mA,  with  conditions  according  to  Fig.  18 
caption. 


216 


T.  Bjazic  et  al.  /  Journal  of  Power  Sources  198  (2012)  203-21 7 


Table  2 

The  parameters  of  the  linear  model  (87)  of  the  boost  converter  supplied  by  the  fuel 
cell  emulator  K,  Tb,  Ti  and  T2  obtained  in  several  operating  points  determined  by  the 
output  current  Iout. 


lout  (A) 

Vin  (V) 

I< 

Tb  (ms) 

h  (ms) 

T2  (ms) 

1 

28.0 

2088 

18.6 

16.0 

74.5 

3 

24.7 

1556 

33.8 

31.7 

22.7 

5 

23.5 

734.6 

94.7 

84.1 

19.5 

7 

22.2 

595.0 

122.4 

102.0 

18.3 

9 

20.5 

461.0 

157.6 

114.2 

17.1 

Fig.  21.  The  Bode  plot  of  the  frequency  characteristics  of  the  obtained  model  (87), 
with  parameters  given  in  Table  2. 

The  transient  responses  in  other  operating  points,  given  in 
Table  2,  are  similar  to  the  responses  in  the  described  boundary 
operating  points,  which  means  that  maximum  error  of  the  linear 
model  is  also  determined  by  the  standard  deviation  of  the  out¬ 
put  voltage  ripple  and  noise.  In  explanation,  the  transfer  functions 
(87),  along  with  the  appropriate  parameters  from  Table  2,  accu¬ 
rately  describe  the  behavior  of  the  peak  current  mode  DC/DC  boost 
converter  in  the  given  operating  point  determined  by  the  output 
current. 

6.  Conclusion 

The  linear  model  of  the  constant  frequency  peak  current  mode 
controlled  DC/DC  boost  converter  supplied  by  the  PEM  fuel  cell 
stack  is  derived  in  this  paper.  The  parameters  of  the  linear  model 
vary  with  the  operating  point,  determined  by  the  load  resistance 
or  output  current.  The  contribution  of  the  paper  is  in  deriving  the 
transfer  function  of  boost  converter  output  voltage  with  respect 
to  the  control  (referent)  current  of  the  converter,  with  included 
compensation  ramp  in  the  continuous  and  especially  in  the  discon¬ 
tinuous  conduction  mode.  Even  though  it  is  not  necessary  for  the 
stabilization  of  the  system  in  the  discontinuous  conduction  mode, 
it  is  hard  to  turn  off  the  compensation  ramp  when  transferring 
from  the  continuous  to  the  discontinuous  conduction  mode  and 
vise  versa,  so  it  must  be  included  in  the  discontinuous  conduction 
mode  as  well.  The  model  is  valid  in  the  whole  range  of  converter 
and  fuel  cell  normal  operating  conditions,  which  makes  it  suitable 
for  advanced  controller  design,  which  will  result  in  approximately 
the  same  dynamic  behavior  of  the  system  irrespective  of  the  oper¬ 
ating  point  and/or  conduction  mode.  The  parameters  of  the  model 
are  obtained  experimentally  in  two  boundary  operating  points  that 
represent  minimal  and  nominal  load  power  consumption,  and  in 
three  other  between  the  boundary  operating  points.  The  linear 


model  errors  are  determined  by  the  standard  deviation  of  the  out¬ 
put  voltage  ripple  and  noise,  so  the  derived  model  with  appropriate 
parameters  accurately  describes  the  system  in  the  given  operating 
point. 
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